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Two fundamental challenging problems of laboratory and astrophysical plasmas are the 
understanding of the relaxation of a cohisionless plasmas with nearly isotropic velocity distribution 
functions and the resultant state of nearly equipartition energy density with electromagnetic plasma 
turbulence. Here, we present the results of a study which shows the role that higher-order-modes 
play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. 

Our main results show that for a thermal plasma the magnetic fluctuations are conflned by regions 
that are bounded by the least-damped higher order modes. We further show that the zone where 
the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer 
wavelengths as the Pe increases. This merging zone has been interpreted as the beginning of the 
region where the whistler-cyclotron waves losses their identity and become heavily damped while 
merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, 
the higher-order modes do not conflne the fluctuations due to the effective higher-temperature 
effects and the excess of suprathermal plasma particles. The analysis presented here considers the 
second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, 
with wave vectors parallel to the uniform background magnetic held, in a flnite temperature 
isotropic bi-Maxwehian and Tsahis-kappa-hke magnetized electron-proton plasma. Our results 
indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these 
quasi modes suggesting that such modes play an important role in the emission and absorption of 
electromagnetic fluctuations in thermal or quasi-thermal plasmas. ©2014 AIP Publishing LLC. 
[http://dx.doi.Org/10.1063/l.4861865] 


I. INTRODUCTION 

The processes of kinetic plasma relaxation and turbu- 
lence generation are responsible for many of the observed 
properties of laboratory and astrophysical plasmas. In partic- 
ular, these processes play an important role in the solar wind 
plasma which is the only cohisionless plasma accessible to 
detailed in-situ satellite observations. The morphology of 
the plasma velocity distribution function (VDF) is generally 
determined by the global structure of the magnetic and 
electric flelds, the nature of the scattering processes (e.g.. 
Coulomb collisions, wave-particle interactions) between 
the source and the measurement location and the boundary 
conditions from which the particles emanate. Detailed meas- 
urements of the VDFs can then provide information about 
the propagation conditions and to some extent about the con- 
ditions in the source region itself. 

In spite of the absence of free-energy for plasma 
instabilities, a flnite-temperature plasma sustains a small but 
detectable spontaneous fluctuations level of absorption and 
emission of various collective modes of oscillations, which 
arises from the discreteness of the plasma particles, in 
other words, from the charge and current fluctuations in 
the plasma. In linear and spatially homogeneous plasmas. 


these spontaneous emissions"^“^ are intimately linked to 
electric and magnetic held fluctuations by the well-known 
fluctuation-dissipation theorem. Spectral properties of the 
scattered fluctuating emission can provide substantial infor- 
mation about the state of a plasma. The analysis of electro- 
magnetic wave scattering is, in fact, one of the most efficient 
methods of plasma diagnostics used in laboratory fusion 
research devices and in space-plasma measurement (e.g.. 
Refs. 13-15). 

The spontaneous thermal emissions of a particular 
eigenmode of a plasma serve as a seed perturbation to drive 
plasma instabilities of the corresponding wave modes when- 
ever some free-energy process becomes available. The study 
of these fluctuations leads also to the concept of correlation 
functions (particularly in time) which play an important role 
in relating the dissipative properties of a system, such as 
the viscous resistance of a fluid or the electrical resistant of 
a conductor, with the microscopic properties of the system 
in a state of near-equilibrium.^^ A good example of the 
fluctuation-dissipation theorem is the quasi-thermal noise in 
space plasmas which is directly linked to fluctuations of 
charged-particle densities and electromagnetic flelds, giving 
rise to scattered radiation. For systems in nearly thermal 
equilibrium, the emitted fluctuating radiation takes on the 
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spectral form of black-body radiation. Because of fluctua- 
tions of charge, current, and electromagnetic fields, the 
particles can be excited, for example, by Langmuir, whistler- 
cyclotron (wc), or Alfven waves, leading to a system that 
cannot be in thermal equilibrium, so that the emitted radia- 
tion deviates from the black-body approximation. However, 
in most cases, the state of the system can be described by the 
linear response to the perturbations from the equilibrium 
state via the fluctuation-dissipation theorem. 

The most striking feature of this theorem is in a funda- 
mental manner, which it relates the fluctuations of a physical 
quantity pertaining to the equilibrium state of a given system 
to a dissipative irreversible process, which in practice is 
realized only when the system is subject to an external per- 
turbation that drives it away from its state of equilibrium. 
Consequently, this theorem enables us to determine the near- 
equilibrium properties of a given statistical system on the 
basis of the knowledge of the thermal fluctuations occurring 
in the system when the system is in one of its equilibrium 
states. Thus, the statistical mechanics of irreversible proc- 
esses is, in a sense, reduced to the statistical mechanics 
of near-equilibrium states (provided we are capable of han- 
dling the time dependent fluctuation processes occurring in 
the system). The fundamental theory of fluctuations was 
advanced by Callen and Welton^^ and later improved and 
refined by Sitenko^^ which provided a detailed account of 
the theory. Although the fundamentals ideas of the theory of 
electromagnetic fluctuations have already been developed 
in these references (Callen and Welton^^ and Sitenko^^), 
the actual application of the formulation to the problem 
of spontaneous emission of whistler-cyclotron waves and 
its connection to the generally ignored higher-order modes 
(horn) have received little attention. 

The electromagnetic dispersion relation of whistler- 
cyclotron waves and ubiquitous higher-order modes in a ther- 
mal plasma has been previously investigated by Gitomer^^ 
and Gitomer and Forslund^^ in a study of collisionless cyclo- 
tron damping of these wave modes, using a numerical linear 
dispersion analysis and non-linear particle simulations in an 
attempt to explain electron heating via these wave eigenmo- 
des. Murata et also investigated similar higher-order 
modes but associated with electrostatic ion-cyclotron waves. 
However, it has not been until recently that these higher- 
order modes have received more attention, in particular, with 
regards to their role in the spontaneous emission of electro- 
magnetic fluctuations in a plasma.^ 

In this paper, we examine the kinetic theory of electro- 
magnetic fluctuations associated with the electron whistler- 
cyclotron waves in plasmas following the theoretical 
approach of Sitenko.^^ In Sec. II, we cover the linear disper- 
sion analysis for parallel fluctuating whistler-cyclotron 
waves including the dispersion characteristics of the higher- 
order modes in a stable-equilibrium plasma. This section 
gives a brief description of the complete set of whistler wave 
modes at parallel propagation, with particular attention paid 
to the topology of the often ignored higher-order modes on 
the dispersion analysis. In Sec. Ill, we discuss the general 
second-order nonlinear theory of fluctuations focusing our 
attention on the emission or absorption of electromagnetic 


fluctuations and their relation to the higher-order modes. The 
general discussion is carried out for two model distribution 
function: (a) bi-Maxwellian and (b) Tsallis-kappa-like 
VDFs.^^“^^ In Sec. IV, we complemented the analytical and 
numerical approaches with a full particle simulation study in 
the stable regime of a plasma under different thermal condi- 
tions. Finally, a discussion of the results and the implication 
to the limiting role of the higher-order modes on the fluctua- 
tions of solar wind plasma observations is presented in 
Sec. V. 


II. LINEAR KINETIC THEORY AND HIGHER-ORDER 
MODES 

Whistler-cyclotron waves are fundamental transverse 
electromagnetic normal modes of magnetized plasmas with 
frequencies between the ion and electron gyrofrequencies, 
CO I I . They are characterized by a helical right- 
handed polarization, which makes them prone to resonate 
with energetic electrons via cyclotron resonance. For that 
reason, they have traditionally received much attention in 
space plasma physics. These wave modes can be driven by 
charged particle thermal anisotropies or by differential 
streaming among various electron populations. For a proper 
understanding of the linear and nonlinear instabilities involv- 
ing such waves, it is essential to understand their mode struc- 
ture (both normal as well as their higher-order modes) 
to develop the ability to recognize them in laboratory and 
space observations. Higher-order modes have been studied 
for electrostatic waves by Derfier,^"^ Murata,^^ and Gary and 
Tokar,^^ and for electromagnetic waves by Gitomer, 
Matsuda,^^ and Astudillo.^^ 

In this section, we describe the linear Vlasov-Maxwell 
theory to determine the full set of electron wave modes 
and higher-order modes for a homogeneous and magnetized 
plasma. We assume electromagnetic fiuctuations propagating 
either parallel or anti-parallel to the ambient magnetic field 
Bq. The linearized dispersion relation tensor Aij{co,k) is 
obtained from the Fourier-Laplace transformed Vlasov- 
Maxwell equations yielding 

Ay(®, k) = Sij{0), k) - ^ • (1) 


The electromagnetic plasma medium properties are 
completely described in terms of the plasma dielectric per- 
mittivity tensor which depends on the complex frequency 
CO and the wave- vector k. The dispersion relation that deter- 
mines the frequencies co = co(k) is given from the condition 
that det(A) = 0. For propagation parallel to Bq, we obtain 
det(A) = A+A_A|| = 0, so that the dispersion relation for 
longitudinal and transverse waves are Ay = 0 and A+ = 0, 
respectively. 

The transverse contributions are given by 


A±(®,fc||) = 1 2^ 

co^ 



.3 

d V ^ 

CO - k\\V\\ ± Qs 


( 2 ) 


and the longitudinal contribution by 
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A||(co,^ll) — 1 + 


^ CO . 


d^v 


V\\ OF, 
CO — k\\v\\ dv\\ ’ 
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where we have defined 


= (co 


dFs 

^11*^11)^ + ^11 *^-L 


dv\\ 


(4) 


Here, co^^ = Anrisq] / rris is the squared plasma frequency and 
= qsBo/(msC^ is the cyclotron frequency, with qs, and 
Ms the density, charge, and mass of the species respec- 
tively, and c is the speed of light. 

For a plasma composed of parallel drifting electrons and 
protons that obey the bi-Maxwellian particle velocity distri- 
bution function 




(i’ll 


- £/ii. 




(5) 


the resultant transverse electromagnetic dispersion relation 
becomes 
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infinite number of root-solutions that are associated to 
higher-order modes. These higher-order modes are also solu- 
tions of the dispersion relation. In the past, these roots have 
been considered unimportant because they are heavily 
damped. As a result, there exists a common interpretation 
that such heavily damped modes do not contribute signifi- 
cantly to the propagation characteristics of the plasma or that 
these modes are not physically realizable in a physical sys- 
tem. In this paper, however, we show that this interpretation 
is misleading. To investigate the nature of these higher-order 
modes, we numerically solved the transverse dispersion rela- 
tion given in Eq. (6). We consider here the role played by 
these higher-order modes and their manifestation in an 
electron-proton plasma. 

The solutions of the transverse dispersion relation 
A+ = 0 are obtained numerically for various plasma condi- 
tions and for the frequency range of interest (i.e., |Qp| < co 
< |Qg|) by considering a non-drifting (^i.e., t/||^ = 0) iso- 
tropic plasma. For the frequency range of interest, the 
dynamics of the ions plays little role since they are more 
massive than electrons and as such high frequencies they 
can almost be considered unmagnetized. In Fig. 1, we 
show the numerical results of the transverse dispersion 
relation (A+ =0) for a warm plasma composed of isotropic 
bi-Maxwellian ions and electrons (T^s/^||s = 1)^ with 
Pe — Pp — ^pet^ie = 1. Here, we show all four- 

quadrants for the wave vector range of —3 < k\\)ie < ?> and 
the normalized frequency — 3 < Re (m)/|Qe| < 3 , where 
ke = c/cope is the electron inertial length but only for the 
positive helicity case (i.e., the positive sign selection in 
Eq. (6). In the top panel [Fig. 1(a)], the solid curves crossing 


and the resultant longitudinal electrostatic dispersion relation 
gives 

A|l(«,^ll) = l-^^ZU), (7) 

s=p,e^\\^\\s 

where is a measure of the temperature anisotropy, 
a||^ = ^J2KBT\\Jms and = \/2kbT\_sI^s are the parallel 
and perpendicular thermal velocities, respectively. Z(c^) and 
Z'((^) are the traditional Fried-Conte plasma dispersion func- 
tion and its derivative with respect to (^, respectively.^^ 

It is well known that linear instabilities can be studied in 
great detail from the solutions of the electromagnetic kinetic 
dispersion relation without any restriction on wavelength, 
growth rate, direction of propagation, etc. Many characteris- 
tics of a given mode (e.g., growth or damping rate, polariza- 
tion, wave vector, compressibility, etc.) can be determined 
by solving the dispersion equation (analytically or numeri- 
cally). However, linear dispersion relation studies tradition- 
ally do not fully take into consideration the topology of 
the whole co-k space. In general, only the fastest growing 
wave and the least damped wave modes are studied by 
following the solutions of the dispersion relation in the 
complex-m plane for a set of fixed macroscopic parameters. 
Nonetheless, any kinetic dispersion relation that includes 
the Fried-Conte plasma dispersion function^^ contains an 



Cfc||/Wpe 
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FIG. 1. (a) Real part of the solution of the transverse dispersion relation A+ 
= 0 for an isotropic Maxwelhan electron-proton plasma, with = 1836, 

cOpej^e = and = 0.01. Solid curves correspond to whistler (wc), ordi- 
nary (O) and extraordinary (X) normal modes, and dotted lines to horn, 
(b) Corresponding imaginary part of the wave modes shown in panel (a). 
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the origin of the (Re ( co ), ^) -plane represent the standard 
forward (third quadrants) and backward (fourth quadrants) 
propagating right-handed wc modes which approaches reso- 
nance at the electron cyclotron frequency -|Qe|. There are 
two other wave modes corresponding to the left-handed 
ordinary (O) and the right-handed extraordinary (X) wave 
modes. Both modes approach asymptotically to the speed 
of light at higher ^|| -values. The ordinary left-handed 
wave modes are the dashed curves that hit the cut-off 
frequency at ^||^ = 0 just below the electron cyclotron 
frequency (m^ + |Qe|), whereas the extraordinary right- 
handed waves are the dashed black curves that hit the cut-off 
frequency at = ^ above the electron cyclotron fre- 
quency (m^ — |Qe|). These waves have their corresponding 
forward and backward propagation counterparts as shown on 
the other quadrants. Fig. 1(a) also includes eight more solu- 
tions of the dispersion relation, of an infinite set of mode 
solutions, corresponding to the least damped horn crossing at 
Re(m/Qe)= —1 at = 0 [red solid lines in Fig. 1(a)]. 
These modes has been studied before and they corresponds 
to the (upper) and (lower) higher-order modes of 
Gitomer and Forslund^^ and Murata et Fig. 1(b) shows 
the corresponding damping rate of the main branches of the 
whistler-cyclotron solutions and also of the higher-order 
modes. Recently, Felten et al^ have investigated a kind of 
aperiodic quasi-modes in an unmagnetized relativistic 
plasma that appear to have similarities with the higher-order 
modes cited here. However, the condition Re{cD) 
required by these aperiodic quasi-modes is not always true 
when compared with the higher-order modes, since they also 
depend on k\\. Nonetheless, some of the higher-order modes 
do satisfy that condition and they could be the aperiodic 
quasi-modes referred by Felten et al^ There have been 
further attempts to solve the dispersion relation to obtain 
approximate analytical solutions to the normal modes for the 
case of low-frequency electromagnetic Alfven waves. 
Such approximations for the higher-order modes have not 
yet been attempted, and here we limit ourselves to the nu- 
merical solution of all these modes. 

Fig. 1(b) shows that for the ordinary and extraordinary 
modes, the damping rate is essentially zero, whereas for the 
whistler-cyclotron branch the damping is nearly zero at long 
wavelengths and then increases substantially as the wave 
vector increases. The damping of the higher-order modes 
increases in the directions given by the red solid lines in 
Fig. 1(b). These modes fill the (Re ( m ), ^) -plane so that the 
least damped higher-order modes form triangle-like zones in 
the dispersion diagram [Fig. 1(a)], with corners at the fre- 
quencies near 0 and — 2Q^ at k\\)ie = ±3. These zones also 
contain portions of the whistler branches dispersion (in the 
fourth and second quadrant). As a consequence, we could 
expect that if free energy is available for the excitation of 
whistler-cyclotron waves, they should emerge in those 
regions bounded by the higher-order modes. In fact, linear 
Vlasov analysis predicts that the whistler-cyclotron branch 
can be unstable with an anisotropic bi-Maxwellian electron 
distribution. As an example, for a temperature ratio 
T^elT\\e = 6 and the same as in Fig. 1 = 0.01), the 

growth rate is positive in the normalized wave vector range 


1.3 < k\\ke < 2.4, reaching its maximum at k\\ke ~ ±1.77 
with Re (cOmax/l^el) ~ —0.03. Although not illustrated here, 
the increase of the temperature anisotropy makes no qualita- 
tive changes in the real part of the dispersion diagram. 

Similar wave mode calculations were carried out for 
larger values of Pe. Figs. 2 and 3 illustrate the solutions of 
A+ = 0 for P^ = O.l and P^ = 0.3, respectively. Note that 
the zones of higher-order damped modes broaden (see Figs. 
2(a) and 3(a)) with increasing P^, whereas the damping rate 
zones decreases (see Figs. 2(b) and 3(b)). Simultaneously, 
the intersection between the least damped mode and the 
whistler-cyclotron mode shifts towards smaller values of 
k\\ke- This result strongly suggests that with increasing P^, 
whistler wave-particle scattering will progressively tend to 
occupy the lower part of the /^-spectrum. As Tu/T\\e 
increases, the whistler branch becomes unstable (not shown 
here) with the positive growth region shifted to smaller wave 
vectors, but still constrained by the least damped 
higher-order modes at small values of k\\kg, and cyclotron 
damping at larger wave vectors. As in the previous case, the 
dispersion diagrams show no qualitative differences between 
the electron isotropy and moderate unstable cases 
(T±e/T\\e < 3). In Sec. Ill, we will discuss an additional ki- 
netic effect that is closely related to the behavior of the 
higher-order modes and the triangle-like zones. The solu- 
tions presented in Figs. 1 to 3 have a mirror image counter- 
parts, when the negative helicity (i.e., negative sign) is used 
in the dispersion relation of Eq. (6). 

III. ELECTROMAGNETIC SPONTANEOUS 
FLUCTUATIONS IN A PLASMA 

The knowledge of the fluctuation spectrum is the start- 
ing point for the computation of the emission of transverse 
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FIG. 2. Similar as in Fig. 1 but with = O.l. 
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FIG. 3. Same as in Fig. 1 but with =0.3. 


electromagnetic waves in a stable thermal equilibrium 
plasma. Our aim in this section is to compute the spectrum 
of these fluctuations, first when the plasma is represented 
by a Maxwellian thermal distributions and later for a more 
general case of a nonthermal distribution modeled by a 
Tsallis-kappa-like distribution function. To describe the 
state of the plasma, we shall use the statistical ensemble 
average of the magnetic field through their deviations from 
the mean values (i.e., fluctuations) which are governed by 
the state of the plasma. The mean value of the fluctuations 
are equal to zero by definition; so the fluctuations are 
described by means of correlation functions (i.e., non-linear 
second order correlators) in the m-k Fourier domain. These 
are average products of fluctuations of different wavelengths 
at different frequencies. The averaging procedure is carried 
out over the mechanical or electrical state of the system and 
over the statistical distribution of various states represented 
by the correlation function. If the plasma medium is 
spatially homogeneous and only stationary states of the 
system are involved in the correlation, then the quadratic 
wavevector-frequency correlation function depends only on 
the relative distance and on the absolute value of the time 
segment between the points at which the fluctuations are 
examined. 

Following standard second-order methods, the classical 
spectral distribution of electric held fluctuations in a magne- 
tized isothermal, and homogeneous plasma can be generally 
expressed (see Sitenko^^) in terms of the fluctuation- 
dissipation theorem where are the components of the 
inverse of the dispersion tensor [Eq. (1)], and ( • • •) repre- 
sents an ensemble average 




( 8 ) 


For parallel propagation and using the Maxwell-Faraday 
equation, we can express this equation in terms of the mag- 
netic fluctuations as 

(9) 

whereas (|^^|| |^) =0, and A+ is given by Eq. (2) above. 

Equation (9) has been solved numerically for the spec- 
trum of magnetic fluctuations as a function of wave vector 
and frequency. The results of this calculation are illustrated 
in Eigs. 4-6 for different j^^-values. Eigs. 4-6 are to be com- 
pared with the dispersion counterparts presented earlier in 
Eigs. 1-3. Eig. 4 displays the contours and intensity scale of 
the magnetic fluctuations obtained from Eq. (9) for the same 
parameters as in Eig. 1. Also shown are the normal mode 
solutions calculated from the dispersion relation A+ = 0, 
representing the whistler-cyclotron, the extraordinary, and 
the ordinary wave normal modes of the system. The color 
scale shows the normalized power levels in logarithmic scale 
between 10“^ and 10“^. Comparing Eigs. 4 and 1, it is clear 
that the spontaneous thermal fluctuations emerge in the 
triangle-like zones, i.e., they are constrained by the least 
damped higher-order modes (dashed lines), which are also 
solutions of the dispersion relation, indicating the zones 
bounding (limiting) the fluctuations of the whistler-cyclotron 
branch. Within these zones, the whistler-cyclotron, although 
still distinguishable from the background fluctuations, is 
heavily damped beyond ^ 1-5 as indicated in Eig. 1(b). 
In addition, Eigs. 5 and 6 show the magnetic fluctuations for 
the cases with =0.1 and = 0.3, respectively. As Pe 
increases and the damped zones become broader, the mag- 
netic fluctuations All the extended regions at frequencies 
even larger than the cyclotron frequency, but otherwise with 
decreasing power. Most of the power is localized along the 
dispersion branch of the whistler-cyclotron mode, where the 
collective properties of the plasma reinforce the spontaneous 
fluctuations, i.e., the individual behavior of electrons. The 
power around the other two eigenmodes, ordinary and extra- 
ordinary waves, is also enhanced but at a signiflcantly lower 
lever than for the whistler-cyclotron waves. This could be 


-1 

-2 

-3 


-4 



FIG. 4. Contours plot and intensity scale of the magnetic fluctuations 
obtained from Eq. (9) for a bi-Maxwellian plasma with the similar parame- 
ters as in Fig. 1. Dotted lines superposed are the least damped higher-order 
modes that conflne the magnetic fluctuations. 
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due to the fact that we use a non-relativistic treatment of the 
fluctuations, and it could be possible that in the relativistic 
case, the O and X modes, which are mostly fast waves (phase 
speed higher than the speed of light in vacuo), should be 
able to resonate with the electrons. Note that, except for the 
symmetry breaking effects of temperature anisotropy or dif- 
ferential streaming due to collective effects, the fluctuation 
power is symmetric around the cyclotron frequency — |Qe| in 
a stable thermal equilibrium plasma. Therefore, for a 
Maxwellian plasma, the least damped higher-order modes 
provides a confining range of whistler-cyclotron fiuctuations 
of the system. These results also shows that as increases 
the point where the whistler-cyclotron waves intersect 
the electromagnetic fiuctuations shift towards longer wave- 
lengths (i.e., smaller ^||^). This intersection region is inter- 
preted as the location where the whistler-cyclotron waves 
begin to loose their identity by merging with the fiuctuations 
in the plasma, which are heavily damped. 


A. Suprathermal effects 

Here, we consider the topology of the spectrum of 
magnetic fiuctuations for a more general case of a nonther- 
mal distribution modeled by a Tsallis-kappa-like distribu- 
tion function.^^“^^ It is assumed that the plasma is 
composed of parallel drifting electrons and protons that 
obey the Tsallis-kappa-like particle velocity distribution 
function given by 







■+i 

K \ di 



1 -(^+1) 


1 + l) 

k'^/^T{k - 1 / 2 ) ’ 


( 10 ) 


Here, the parameter k represents the degree on non- thermal 
effects. When k is small, the velocity distribution in Eq. (10) 
represents an effectively hotter distribution than a 
Maxwellian, whereas for k large and approaching infinity, 
the effective temperature reduces to that of a Maxwellian 
distribution. This effective temperature of a Tsallis- 
kappa-like distribution can be defined in terms of the 
Maxwellian temperature T by obtaining the second moment 
of the distribution function Eq. (10) to yield 


T, 


K 

k-312 


T. 


(11) 


Substituting Eq. (10) into Eqs. (2)-(4), we can write the 
perpendicular dispersion tensor as follows: 


klc'^ 

s=p,e 


= 


r(K) 


7rl/2K'/2r(K- 1/2) 

+“ dt r, f 
1 +- 
K. 


-oo (t - 0 


, lm(^)>0, (12) 


where, for k> 1/2, the function Zk{^) can be written in 
terms of the Gauss hypergeometric function 2 F 1 


' ^3/2 2^1 


1, 2k; K + 1; 1 ( 1 - 


(13) 


and the arguments and have been defined in Eq. (6) 
above. Here, Z^((^) is the analogous plasma dispersion func- 
tion for a kappa-like distribution as defined by Mace and 
Hellberg^^ and Mace and Sydora.^^ 

We now present an estimate of the power spectrum of 
magnetic fiuctuations as defined by Sitenko,^^ but solving 
the dispersion relation for the Tsallis-kappa-like velocity 
distribution in Eq. (10), using Eqs. (9) and (12). Eig. 7 repre- 
sents the numerical results of the transverse dispersion 
relation (A+ = 0) for a warm plasma composed of a non- 
drifting (^i.e., t/||^ = 0) isotropic Tsallis-kappa-like velocity 
distribution of ions and electrons = 1), with 

Pe = Pp = 0.1 and (J0pe/\kle\ = 1, which are similar parame- 
ters as those used in Eig. 2 for a bi-Maxwellian plasma. Here 
we use a value of k: = 3 that represents a particle velocity 
distribution with a substantial suprathermal tail. For compar- 
ison purposes, here, we show the same range of wave vectors 
3 3 and normalized frequency 3 R.e {«)/ 

|Qe| < 3 shown in Fig. 2. The solid curves of Fig. 7(a) cross- 
ing the origin of the (Re (m), /:) -plane represent the standard 
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Cfcll/Wpe 

FIG. 7. (a) Real part of the solution of the transverse dispersion relation 
A+ — 0 for an isotropic electron-proton plasma following nonthermal 
distribution modeled by a Tsallis-kappa-like distribution function with 
k: = 3, rriplme = 1836, (OpelO^e = 1, and = 0.1. Solid curves correspond 
to whistler (wc), ordinary (O) and extraordinary (X) normal modes, and 
dotted lines are the horn, (b) Corresponding imaginary parts of wave modes 
shown in panel (a). 

forward (third quadrant) and backward (fourth quadrant) 
propagating right-handed whistler-cyclotron modes which 
approaches resonance at the electron cyclotron frequency 
— |Qg|. The other two wave modes correspond to the left- 
handed ordinary and the right-handed extraordinary wave 
modes. Both modes approach asymptotically the speed of 
light at higher /^-values. In Fig. 7(a), the ordinary left-handed 
wave modes are the black-dashed curves that hits the cut-off 
frequency at = 0 just below the electron cyclotron 
frequency (m^|Qe|), whereas the extraordinary right- 
handed waves are the black-dashed curves that hits the cut- 
off frequency dX k\\ke above the electron cyclotron fre- 
quency (m^ — |Qe|), similar to Fig. 2(a) above. These waves 
have their corresponding forward and backward propagation 
counterparts as shown on the other quadrants. Fig. 7(a) also 
includes two more solutions of the dispersion relation, corre- 
sponding to the least damped higher-order modes crossing at 
Re(m/Qe)=— 1 at = 0 [red solid lines in Fig. 7(a)]. ^ 

Although Figs. 7 and 2 have many similarities, one striking — 

feature between the two solutions is that the number of 3 
higher-order modes is finite for a Tsallis-kappa-like distribu- 
tion. Furthermore, the total number of higher-order modes 
n obeys the relationship n = {k — l). Thus, for the case 
depicted in Fig. 7 where k = 3, the total number of higher- 
order modes is four. The higher-order modes in Fig. 7(a) are 
grouped very close together due to the rather effective higher 
temperature of the non-thermal distribution compared to 
the Maxwellian case. The Z^-function used in our calcula- 
tions is essentially the same one defined by Hellberg and 


Mace^^ in their Eq. (10). Using a similar treatment as 
Summers and Thorne,^^ it can be shown that the dispersion 
relation is a polynomial of order k: + 2 when k is integer and 
thus admits k: + 2 complex solutions. After removing the 
electromagnetic modes O, X, and whistler, then the net num- 
ber of higher-order modes reduces to k: — 1 . Alternatively, it 
can be shown that the 2^1 hypergeometric function for k in- 
teger is represented by a finite-order polynomial, whereas for 
K real, it has branch-cuts with infinite solutions. 

Fig. 7(b) shows the damping rate of the main branches 
of the normal mode solutions and of the higher-order modes. 
For the whistler-cyclotron branch, the damping is nearly 
zero at long wavelengths and then increases substantially as 
the wave vector increases. For the ordinary branch, the 
damping is greater that for the Maxwellian case. This is 
because for the Tsallis-kappa-like distribution the effective 
temperature for k = 3 is higher than the Maxwellian, and 
therefore such cases are quite similar to a higher temperature 
Maxwellian case. For the extraordinary modes, the damping 
rate remains essentially zero. The damping of the higher- 
order modes increases in the directions given by the red- 
solid lines in Fig. 7(b) but the triangle-like zones formed 
by these modes in the (Im(m), ^)-plane show smaller disper- 
sion among themselves than in the Maxwellian situation for 
the same = 0. 1 . However, they seem to be similar to the 
higher- case of the Maxwellian distribution [depicted in 
Fig. 3(b)]. 

In addition. Fig. 8 shows the corresponding electromag- 
netic fluctuations for the case of = 0.1 of a Tsallis-kappa- 
like velocity distribution function with k = 3. As before, it 
shows clearly the presence of the whistler-cyclotron, extraor- 
dinary and ordinary normal modes of the system. The color 
scale shows logarithmic normalized power levels between 
10“^ and 10“^. Note that in this case, the electromagnetic fluc- 
tuations are not bounded anymore by the least-damped higher 
order modes. This is because the effective temperature of this 
plasma [given by Eq. (11)] is effectively higher than a 
Maxwellian, i.e., 7^ = 27 for such ^-value. This is expected 
since the Tsallis-kappa-like distribution is a non-thermal dis- 
tribution with an excess of suprathermal particles providing 
such effective higher temperatures. Thus, the electromagnetic 
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FIG. 8. Contours plots and intensity scale of the magnetic fluctuations 
obtained from Eq. (9) for a Tsallis-kappa-like plasma (k — 3) with the same 
parameters as in Fig. 7 including the least damped higher-order modes 
(dashed lines). 
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fluctuations can grow outside the limiting region of the least 
damped higher-order modes. Although not shown in here, 
this is also valid for different values of and k. Thus, in 
summary, the higher order modes do not limit the electro- 
magnetic fluctuations in a non-thermal plasma. This signature 
can be used as a diagnostic for detecting the presence 
of non-thermal plasmas in plasma systems where we have 
access to magnetic fluctuation data but we do not have direct 
in-situ measurements of the distribution function. 

IV. FULL PARTICLE SIMULATIONS 

We now compare our theoretical calculations against 
full particle simulations for the same set of parameters used 
in Secs. I-III. We carry out particle-in-cell (PIC) simulations 
where the ions and electrons are treated kinetically and self 
self-consistently with the electric and magnetic field fluctua- 
tions. The simulation code is analogous to the KEMPOl 
code of Omura and Matsumoto^^ except for the schemes to 
correct the longitudinal electric field, the interpolation to the 
particle position and velocities, and the normalization. The 
simulation code is initiated with a Maxwellian distribution 
with a density and temperature consistent to the parameters 
used in Figs. 1-3. Periodic boundary conditions are imposed. 
The equations for the particle position and velocities are 
solved using the Boris-Buneman scheme and evolved in 
time by means of the leapfrog method. Moment charge 
densities and currents are calculated by a second-order 
Gaussian-like shape function by interpolation to the particle 
positions and velocities followed by accumulation of their 
charge and current contributions for each cell. This approach 
reduced considerably the noise level to obtain cleaner and 
smoothed results. Maxwell’s equations for the electric and 
magnetic fields are solved via a second-order finite differ- 
ence scheme. The longitudinal electric-field is corrected 
according to the Marder^^ and Langdon^^ schemes to assure 
that the longitudinal electric field determined via the 
Ampere’s equation is self-consistent with that determined 
from Poisson’s equation throughout the simulation. Time is 
normalized in units of the electron cyclotron frequency, 
particle positions are normalized to the electron inertial 
lengths, whereas all velocities are normalized to the speed- 
of-light. Because the dimensional system of equations was 
written in cgs-units, the electric and magnetic fields were 
normalized similarly. 

The temporal collisionless cyclotron damping of whis- 
tler waves has been studied by Scarf, Ossakow et al.^^ 
Gitomer,^^ and Gitomer and Forslund^^ using numerical sim- 
ulations. The calculated numerical dispersion of these modes 
agrees well with the prediction of a finite temperature linear 
dispersion theory. In our 1.5-D (i.e., one spatial and 3 veloc- 
ity dimensions) simulations, the plasma is presumed to be 
collisionless, homogeneous, and magnetized, and it consists 
of a single ion and electron components with realistic mass 
ratio mplrue = 1836. Initially, electrons and protons satisfy 
an isotropic Maxwellian distribution, so that no free energy 
is present and the plasma is stable. To avoid side effects 
such as multi beam formation at low no quiet start 
particle-loading scheme was employed. The number of grid 


cells is 1024 with 1000 particles per species (electrons 
and protons) per cell. Similar results are obtained if 1600 
particles per cell are used. The length of the simulation box 
is 512 in unit of the electron inertial length 2^, the time step 
is Q^Aa = 0.05, and periodic boundary conditions were 
imposed for particles and fields. In this section, we focused 
only on the simulation results generated by a Maxwellian 
velocity distribution. The simulation results for an initial 
Tsallis-kappa like particle distribution will be presented later 
in a subsequent paper. 

Figs. 9-11 represent the simulation counterparts of the 
theoretical results presented earlier in Figs. 4-6. Fig. 9 dis- 
plays the co-k power spectrum of the magnetic field fluctua- 
tions for = 0.01, obtained by Fourier-transforming the 
transverse magnetic field data both in space and time. The 
color scales show the power levels (normalized to the maxi- 
mum wave power) in logarithmic scale between 10“^ and 
10~^. After decomposing the wave fluctuations into positive- 
and negative-helicity parts, we can clearly recognize the for- 
ward (backward) propagating right-hand modes in the third 
(fourth) quadrants, respectively. Fig. 9 also shows details of 
the spontaneously emitted magnetic fluctuations, in agree- 
ment with Figs. 1 and 4. Note the triangle-shape zone in 
Fig. 9 beginning at m/|Qe| = —1 near \k\\ke\ =0 that encom- 
pass the whistler-cyclotron branch beyond hK > 1 . In these 
triangle- shape zones, the whistler-cyclotron branch merges 
with the fluctuations, loosing its identity and becomes heav- 
ily damped. Figs. 10 and 11 depict the situations for the 
cases with = 0.1 and = 0.3, respectively, and illustrate 
that this merging region of the whistler-cyclotron waves and 
the fluctuations shift to lower k\\)ie values as increases. 
Clearly, for the j^^-values investigated, the spontaneous 
magnetic fluctuations are enhanced in the vicinity of the 
whistler-cyclotron modes, which yet quickly lose their 
identity at large values of k\\)ie with increasing Thus, if 
enough free energy is available to overcome the natural 
threshold imposed by the thermal fluctuation level, the 
growth of linear and even nonlinear instabilities could be 
understood in terms of the amplification of some of the fluc- 
tuation modes. A quick comparison of Figs. 9 to 11 with 
Figs. 1 to 3 shows that although we have chosen the positive 
helicity solutions in the simulations, the results also includes 



FIG. 9. Contour plot and intensity scale of the normalized magnetic power 
spectrum obtained from numerical particle-in-cell simulations of a electron- 
proton Maxwellian plasma, with mpjme — 1836, cOpejO^e = and — 0.1 
for an initial bi-Maxwellian plasma. 
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FIG. 10. Similar to Fig. 9 but with = 0.1. 



the mirror-image counter parts, corresponding to negative 
helicity, but at a much lower intensity. This is consistent 
with the mirror-image solutions obtained by changing the 
helicity sign in Figs. 1 to 3. 

Fig. 12 shows a characteristic power spectrum obtained 
from the numerical PIC- simulations. The plot corresponds 
to a spatial Fourier mode = 0.075 for the case of 
= 0.01 and C0pe/\Q^e\ = 1- The principal peak is the 
whistler cyclotron wave with a real frequency at co/\Q.e\ 
= —0.31. The lower amplitude higher frequency modes near 
m/|Q^| = — 1.0 are the fluctuations previously discussed. 
These results are consistent with the result plots obtained by 
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Gitomer and Forslund^^ (see their Fig. 3). There are two 
other peaks associated with the ordinary (O) and extraordi- 
nary (X) wave modes present in our simulation that do not 
appear in Gitomer and Forslung^^ paper due to the fact that 
they use an implicit electromagnetic PIC-code that Alter all 
the speed-of-light wave modes, whereas we have used and 
explicit PIC-code. Nonetheless, our results are consistently 
similar. 

V. CONCLUSIONS 

A description of the role of higher-order damped modes 
and their relation to the spontaneous magnetic fluctuations in 
magnetized thermal and nonthermal plasma has been pre- 
sented. We have shown the relevance of higher-order modes in 
the emergence of spontaneous plasma fluctuations due to 
quasi-thermal noise through the fluctuation-dissipation theo- 
rem. From kinetic theory, the spectral distribution of the mag- 
netic fluctuations is topologically constrained by the structure 
of the heavily damped higher order modes in a thermal plasma. 
These regions becomes broader as the plasma increases and 
they determine where instabilities could emerge if free energy 
is made available to the system. We also have shown that 
the zone where the whistler-cyclotron waves merges with the 
electromagnetic fluctuations shifts to longer wavelengths as 
)5g-increases. The zones where the whistler-cyclotron branch 
merges with the fluctuations are interpreted as regions where 
these waves lose their identity and becomes heavily damped. 
All our linear dispersion calculations and theoretical spectral 
fluctuation spectrum results exhibit remarkably consistent 
agreement with fully self-consistent numerical simulations for 
the same plasma conditions. 

We have also carried out similar calculations for a non- 
thermal plasma modeled by a Tsallis-kappa-like velocity 
distribution function. In this case, the electromagnetic fluctu- 
ations are unconstrained by the higher-order modes because 
the Tsallis-kappa-like distribution has an excess of suprather- 
mal particles that produce an effective temperature hotter 
than a thermal Maxwellian plasma. This result provides a 
critical observable signature, to determine the nature of a 
plasma, via its velocity distribution function, when using 
electromagnetic fluctuations data to diagnose the plasma 
even when direct in-situ measurements of the distribution 
function are not available. This important signature from ki- 
netic theory should have applications to the understanding of 
the heating of electrons in laboratory and astrophysical plas- 
mas. The calculations presented in this paper can be easily 
extended by analogy to the case of the low-frequency electro- 
magnetic Alfven waves in which the ions play the major role. 
In summary, our results indicate that the spontaneously emit- 
ted electromagnetic fluctuations are in fact enhanced over 
these higher-order modes, suggesting that such modes play 
an important role in the emission and absorption of electro- 
magnetic fluctuations in thermal or quasi-thermal plasmas. 
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